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1/^ ' We have applied a new gauge-invariant, noncompact, Monte Carlo method to simulate the (7(1), SU{2), and 

^\ , SU{3) gauge theories on 8"^ and 12'' lattices. The Creutz ratios of the Wilson loops agree with the exact results 

■ for ?7(1) for /3 > 0.5 apart from a renormalization of the charge. The SU{2) and SU{S) Creutz ratios robustly 
display quark confinement at /3 = 0.5 and (3 = 2, respectively. At much weaker coupling, the 5(7(2) and 5'?7(3) 

^ , Creutz ratios agree with perturbation theory after a renormalization of the coupling constant. For 5(7(3) the 

^ ' scaling window is near /3 = 2, and the relation between the string tension a and our lattice QCD parameter Al 

["T I [ is y/a ^ 5Ai. 
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1. INTRODUCTION 

In compact lattice gauge theory, gauge fields 
are represented by group elements rather than 
by fields, and the action is a periodic function 
of a gauge-invariant lattice field strength. The 
periodicity of the action entails spurious vacua. 
The principal advantage of noncompact actions, 
in which gauge fields are represented by fields, is 
that they avoid multiple vacua. 

The first gauge-invariant noncompact simula- 
tions were carried out by Palumbo, Polikarpov, 
and Veselov |jl[. They saw a confinement signal. 
Their action contains five terms, constructed from 
two invariants, and involves (noncompact) auxil- 
iary fields and an adjustable parameter. 

The present paper describes a test of a new 
way H] of performing gauge-invariant noncom- 
pact simulations. Our action, which is similar to 
one term of Palumbo's action, is exactly invari- 
ant under compact gauge transformations, is a 
natural discretization of the classical Yang-Mills 
action, and reduces to Wilson's action when the 
gauge fields are compactified. In this method 
there are fewer auxiliary fields than in Palumbo's 
method, and they are compact group elements 
representing gauge transformations. 
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We have used this method to simulate U{1), 
SU{2), and SU{i) gauge theories on S'' and 12" 
lattices. The Creutz ratios of Wilson loops agree 
with the exact results for U{1) for (i > 0.5 apart 
from a renormalization of the charge. The 5(7(2) 
and 5(7(3) Creutz ratios clearly show quark con- 
finement at (3 = 0.5 and (3 = 2, respectively. 
At much weaker coupling, the SU{2) and 5(7(3) 
Creutz ratios agree with perturbation theory with 
a renormalized coupling constant. For 5(7(3) 
there is a scaling window near (3 = 2, and the 
string tension cr is related to the lattice QCD pa- 
rameter Al by « 5Al. If « 420 MeV, 
then our Al is about 84 MeV, and at (3 = 2 our 
lattice spacing a is about 0.4 fm. 



2. THE METHOD 

For massless fermions, the continuum action 
density is ipi"/fj,dfj,ip. A suitable discretization of 
this quantity is ii(j{n)jij,[ip{n + e^) — ip(n)]/a in 
which n is a four-vector of integers representing 
an arbitrary vertex of the lattice, is a unit 
vector in the /ith direction, and a is the lattice 
spacing. The product of Fermi fields at the same 
point is gauge invariant as it stands. The other 
product of Fermi fields becomes gauge invariant 
if we insert a matrix (n) of gauge fields 
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that transforms appropriately. Under a gauge 
transformation represented by the group elements 
U{n) and U{n + e^), the required response is 

l+iagA'^{n) = U{Ti)[l+iagAf,{n)]U~\n+ef,).{2) 

Under this gauge transformation, the lattice field 
strength 

+ ig [A^{n)Af_,{n + e^) 

-A^(n)A^(n + e^)], (3) 

which reduces to the continuum Yang-Mills field 
strength in the limit a —t 0, transforms as 

^;.(") = U{n)F^,{n)U-\n + + e,). (4) 

The field strength F^i, (n) is antisymmetric in the 
indices /i and but it is not hermitian. To make 
a positive plaquette action density, we use the 
Hilbert-Schmidt norm of Ff^,j{n) 

S ^ ^TT[Fl{n)F^,{n)], (5) 

in which the generators Ta of the gauge group are 
normalized as Tr(raTb) — kSab- Because Ff^„{Ti) 
transforms covariantly (^, this action density is 
exactly invariant under the noncompact gauge 
transformation (^). 

In general the gauge transformation (||) with 
group element U{n) — exp(— iago^Ta) maps the 
matrix of gauge fields Afj_{n) — TaA'^{Ti) outside 
the Lie algebra, apart from terms of lowest (ze- 
roth) order in the lattice spacing a. We use this 
larger space of matrices. We use the action (||) 
in which the field strength (|^) is defined in terms 
of gauge- field matrices (n) that are the images 
under arbitrary gauge transformations 

AJn) =VA"(n)W-'^ + —V (V-^ -W-'^) (6) 

ag 

of matrices of gauge fields defined in the 

usual way, A°^{n) = TaA°'°(n). The group ele- 
ments V and W associated with the gauge field 
Afj,{n) are unrelated to those associated with the 
neighboring gauge fields A^{n + e^), A^{n), and 
A,y{n + e^). 



The quantity I -I- igaA^{n) is not an element 
L^{n) of the gauge group. But if one compact- 
ified the fields by requiring I -I- igaA^in) to be 
an element of the gauge group, then the matrix 
yip (n) of gauge fields would be related to the link 
L^{n) by A^{n) = {Lf^{n) - l)/{iga), and the ac- 
tion (H) defined in terms of the field strength (^) 
would be, mirabile dictu, Wilson's action: 

k - lR.TvL^{n)L^{n + e^)Ll{n + e^)Ll[n) 

3. RESULTS 

We have tested this method by applying it to 
the C/(l), SU{2), and S'C/(3) gauge theories on 8^ 
and 12* lattices. In most of our initial configu- 
rations, the unitary matrices V and W and the 
hermitian gauge fields A" were randomized. For 
thermalization we allowed 50,000 sweeps for C/(f ), 
10,000 for SU{2), and 100,000 for SU{i). Our 
Wilson loops are ensemble averages of ordered 
products of the binomials 1 -I- iagAfj^{n) rather 
than of the exponentials exp[iagA^{n)] around 
the loop. 

For U{1) and for P — 1/g^ > 0.5, our mea- 
sured Creutz ratios |^ of Wilson loops agree 
with the exact ones apart from finite-size effects 
and a renormalization of the charge. For in- 
stance at /3 = 0.5 on the 12* lattice, we found 
X(2,2) = 0.147(1), x(2,3) = 0.103(1), x(2,4) = 
0.090(1), x(3,3) = 0.049(1), x(3,4) = 0.034(1), 
and x(4,4) = 0.020(2). The first three of these 
x's are equal to the exact Creutz ratios for a 
renormalized value of /3r = 0.45; the last three 
are smaller than the exact ratios for f3r — 0.45 
due to finite-size effects by 3%, 8%, and 16%, re- 
spectively. 

But at stronger coupling, the extra terms 
ig [Ai, {n)A^ (n -I- e^) - yl^ in)A^ {n + e^)] in the 
lattice field strength Ff^^{n) eventually do pro- 
duce a confinement signal. For example, at 
(3 ~ 0.375, our measured Creutz ratios on the 
12* lattice are: x(2, 2) = 0.906(5), x(2,3) = 
0.909(21), x(2,4) = 0.85(10), x(3,3) = 0.62(24), 
and x(3,4) = 0.6(16). 

For SU{2) on the 8* lattice at /3 = A/g^ = 0.5, 
we found x(2,2) = 0.835(3), x(2,3) = 0.852(12), 
X(2,4) = 0.865(60), and x(3,3) = 0.94(23) which 
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within the limited statistics clearly exhibit con- 
finement. At [3 — 1, our six Creutz ratios track 
those of tree-level perturbation theory for a renor- 
malized value oi (3r = 1.75. 

For SU{i) at /3 = 6/g2 = 2 on the 12^ lat- 
tice, we found in ten independent runs x(2, 2) = 
0.838(1), x(2,3) = 0.826(3), x(2,4) = 0.828(13), 
X(3,3) = 0.793(42), x(3,4) = 0.47(25), and 
;^(4, 4) = 1.2(86). Within the statistics, these 
results robustly exhibit confinement. At much 
weaker coupling, our ratios agree with perturba- 
tion theory apart from finite-size effects and after 
a renormalization of the coupling constant. 

4. SCALING 

We used an 8^ lattice to study the scaling of the 
lattice spacing a with the coupling constant g for 
SU{'i). The two-loop result for the dependence of 
the string tension aa? upon the inverse coupling 
13 = 6/g2 is 
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(7) 



If we set ^/a K, (5.0 ± 0.4)Al, then our xihjV^ 
fit this formula for 1.9 < /? < 2.1 as shown in the 
figure. A string tension y/a « 420 MeV implies 
that Ai « 84 MeV, which is about 11 times closer 
to the the continuum A.^^ than is the parameter 
Alw ~ 7.9 MeV of Wilson's method. At /? = 2, 
our lattice spacing a is about 0.4 fm. 
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Figure 1. The SU{3) Creutz ratios xihj): the 
scaling predictions for the string tension aa^ , and 
the tree-level curves for x(2, 2) (dots) and x(4, 4) 
(dashes) are plotted against (3. For 1.9 < /? < 2.1, 
some of the symbols of the x's overlap. 
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